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Beam dynamics at the main LEBT of RAON accelerator 

Hyunchang jiif] and Ji-Ho Jang 
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Abstract 

The high-intensity rare-isotope accelerator (RAON) of the Rare Isotope Science Project (RISP) 
in Daejeon, Korea, has been designed to accelerate multiple-charge-state beams. The ion beams, 
which are generated by Electron Cyclotron Resonance Ion Source (ECR-IS), will be transported 
through the main Low Energy Beam Transport (LEBT) system to the Radio Frequency Quadrupole 
(RFQ). While passing the beams through LEBT, we should keep the transverse beam size and 
longitudinal emittance small. Furthermore, the matching of required twiss parameter at the RFQ 
entrance will be performed by using electro-static quadrupoles at the main LEBT matching section 
which is from the multi-harmonic buncher (MHB) to the entrance of RFQ. We will briefly review 
the new aspects of main LEBT lattice and the beam matching at the main LEBT matching section 
will be presented. In addition, the effects of various errors on the beam orbit and the correction of 
distorted orbit will be discussed. 
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FIG. 1: Layout of RAON heavy ion accelerator. 


I. INTRODUCTION 


The Rare Isotope Science Project (RISP) was established in December 2011, in the Insti¬ 
tute of basic science (IBS) to build a heavy ion accelerator, which is named as RAON (Rare 
isotope Accelerator Of Newness) [THlj in 2012, for various science programs. The project of 
RAON heavy ion accelerator is in progress to produce a variety of high-intensity stable ion 
beams and rare isotope beams to be used for a wide range of basic science researches and 
various applications in astrophysics, atomic physics, material science, nuclear science and so 
on. 

In the RAON heavy ion accelerator, the in-flight fragmentation (IF) and Isotope Sepa¬ 
ration On-Line (ISOL) [S] methods will be used to produce various isotopes. In case of a 
uranium beam, the hnal beam energy and power to the various targets will be 200 MeV/u 
and 400 kW, respectively. To achieve this goal, the IF system uses a diver linac Pd which 
consists of a 28 GHz superconducting electron cyclotron resonance ion source (ECR-IS), a 
main low energy beam transport (LEBT) section |8], a 300 keV/u REQ (radio-frequency 
quadrupole) accelerator, a medium energy beam transport (MEBT) section, a 18.5 MeV/u 
low energy superconducting linac (SCLl), a charge stripper section (CSS), and a high en¬ 
ergy supoerconduction linac (SCL2). Unlike this system, the ISOL system uses a 70 MeV 
cyclotron as the driver to deliver 70 kW beam power to the ISOL target. The rare iso¬ 
tope beams created by the ISOL system, are reaccelerated by a post linac accelerator which 
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FIG. 2: Layout of main LEBT. 

consists of a post LEBT, a post RFQ, a post MEET, and a 16.5 MeV/u superconducting 
linac SCL3. The rare isotope beam accelerated by SCL3 will be delivered to the low energy 
experimental hall or to the SCL2 after passing through P2DT (post-to-driver transport) 
line. The schematic layout of RAON accelerator is shown in Fig. 

Recently, the lattice of main LEBT is partly modified to decrease the transverse beam 
size and improve the longitudinal acceptance at the SCLl. The simplified layout is shown 
in Fig. The roles of main LEBT system are to select targeted heavy-ion beams among 
various ion beams generated by ECR-IS and to deliver these seleted beams to the RFQ with 
keeping the transverse beam size and longitudinal emittance small. Among various charge- 
state beams created by ion source, only one or two charge state beams have to be transported 
to RFQ to match the requirements of beamline experiment. This charge selection is carried 
out between two 90-degree bending magnets which provide horizontally dispersive region. 
In the dispersive section, each beam trajectory is changed by the momentum difference 
of each particle which is directly related to the charge of mass ratio. Selected beams are 
longitudinally bunched by a multi-harmonic buncher (MHB) [101 E] to increase the longi¬ 
tudinal transmission at RFQ and SCLl. If two-charge beams are transported, the difference 
of particle velocities has to be equalized by VE before entering to RFQ. In addition, the 
twiss parameters have to be matched to the required values at the RFQ entrance to increase 
beam transmission. Therefore, electrostatic quadrupoles, which are located between MHB 
and RFQ, are used to match beam twiss parameters to the requirement of RFQ. 

In this paper, we will review the upgraded lattice design of main LEBT and discuss the 
result of beam dynamics simulation with a particle tracking code, TRACK [12]. At first, the 
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FIG. 3: Layout of MHB. 
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FIG. 4: Longitudinal axial field map of MHB. 

new layout and field map of MHB will be discussed. Next, the matching of twiss parameter 
at the entrance of RFQ will be described. Finally, the effects of errors induced by magnets 
throughout the main LEBT will be discussed and the correction of those errors will be 
presented with limited number of BPMs and correctors. 


II. BEAM DYNAMICS IN THE MAIN LEBT 

A. Design of MHB 

To transport two charge-state beams from ECR-IS to RFQ, beams have to be bunched 
before entering the RFQ to increase the transmission rate and longitudinal acceptance in the 
downstream of main LEBT. Eor that reason, there are two bunchers, a MHB and a velocity 
equalizer (VE), at the end of main LEBT. Among two bunchers, the MHB is an effective 
system for bunching low-intensity DC ion beams. The MHB in the main LEBT operates 
with three frequencies (40.625, 81.25, and 121.875 MHz) to produce an approximately linear 
sawtooth in the voltage as a function of time to minimize the longitudinal emittance growth. 
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FIG. 5: Longitudinal particle distributions before MHB (left top), after MHB (right top), after 
VE (left down), and before RFQ (right down). 

The fundamental frequency of MHB is 40.625 MHz which is half of the RFQ frequency. On 
the other hand, the VE equalizes the velocity of two beams as a higher energy beam gets 
lower energy gain than a lower energy beam during passing through VE. The distance 
between two bunchers is given by da 

^ , [2^ VQoiQo - 1 ) 

Lmhb,ve - 

where Xrfq is the wavelength of the RFQ, Vq is the accelerating voltage, A is the mass 
number, Qq is the hightest charge state of ions. 

Recently, the design of MHB is modified newly to improve the longitudinal acceptance at 
SCLl linear accelerator. Figure]^ shows the side cross-section of the MHB. The total length 
of MHB is 20 cm and the aperture diameter at middle point is 1 cm. Therefore a maximum 
transverse beam size should be less than 1cm to avoid beam loss at the MHB. In order 
to achieve this requirement, three electrostatic quadrupoles ahead of MHB are used. The 
voltage amplitude of MHB is calculated depending on the particle velocity and distances 
between MHB, VE, and RFQ. The axial field map of MHB is shown in Fig. 

With newly designed MHB, two charge-state uranium beams, 128U33-I- and 128U34-I-, 
are tracked with longitudinally uniform distribution as shown in Fig. After MHB, beams 
have saw-tooth energy modulation as keeping energy difference 1.5 %. This energy difference 
disappears after passing through the VE. Consequently, the rms longitudinal emittance is 
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FIG. 6: Matching of twiss parameters at RFQ entrance with beam envelop in main LEBT. Last 
five electrostatic quadrupoles are used to match the required twiss parameters at RFQ entrance. 
Transverse and longitudinal beam profiles at RFQ entrance are shown in the figure. 

decreased at the SCLl entrance after modifying the design of MHB. The value of rms 
longitudinal emittance was 0.057 degxMev/u for previous MHB and it decreases to 0.045. 
Furthermore, this emittance is much smaller than the result of without MHB and VE, 0.096. 

B. Matching twiss parameter 

One of important issues in the main LEBT is the matching of twiss parameters at the 
RFQ entrance. For higher beam transmission rate, the specific twiss parameters, a beta- 
function 5.67 cm/rad and an alphafunction 0.324, are required at the RFQ entrance. To 
do that, at least 4 electrostatic quadrupoles are needed ahead of RFQ. However, a max¬ 
imum beam size is close to the chamber inner radius after matching with 4 electrostatic 
quadrupoles. For that reason, 5 electrostatic quadrupoles are used between MHB and RFQ 
to fulfill all requirements because more quadrupoles can not be installed by reason of lim¬ 
ited spaces between MHB and RFQ. Twiss parameters are matched by using simulation 
codes ELEGANT [Hj and TRACK. Finding a matching solution with single particle beam 
is performed using ELEGANT, and then the matching is carried out a multi-particle beam 
using a TRACK code. Figure shows the TRACK simulation result with a uranium beam, 
238U33.5+ after finishing twiss parameter matching with 100,000 macro-particles. As a 
result, the matching of twiss parameters is performed with high confidence, that is less than 
0.1 % error at the RFQ entrance. 
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TABLE I: Tolerance of error sources. 


- 

- 

Value 

Unit 

Initial position 

x/y 

0.3/0.5 

cm 

Initial angle 

xp/yp 

5.0/12.0 

mrad 

Solenoid 

misalignment 

2.0 

cm 

Solenoid 

field amplitude 

0.7 

% 

Bending 

misalignment 

3.0 

mm 

Bending 

tilt 

2.5 

mrad 

Bending 

field amplitude 

0.2 

% 

Quadrupole 

misalignment 

0.6 

mm 

Quadrupole 

tilt 

8.0 

mrad 

Quadrupole 

field amplitude 

0.15 

% 


III. ERROR ANALYSIS 

There are many components which can distort beam orbits in the main LEBT like magnet 
misalignment, field error, tilt and so forth. Each error component is listed in Table |I] and 
the tolerance of those errors is also examined to check the sensitivity of each error source. 
Tolerance is calculted with a particle loss during simulations of 100 random seeds for each 
error. Among them, electrostatic quadrupole errors are a main source of beam distortion 
because there are 20 electrostatic quadrupoles in the main LEBT. 

Subsequently, an error study is performed with given errors as listed in Table |TI] for 100 
random seeds. The distortion of orbit trajectories is shown in Figure The orbit distor¬ 
tion is most seriously appeared between MHB and RFQ because there are 5 electrostatic 
quadrupoles with short drift spaces. In the main LEBT, there is no enough space to install 
lots number of correctors and BPMs, therefore proper number and location of correctors and 
BPMs have to be investigated. In addition, the distortion of beam orbit should be corrected 
strictly at the RFQ entrance. Under these limitations, 2 correctors and 2 BPMs are used 
to correct the distorted orbits. Each corrector is assumed to give a horizontal and vertical 
kicks at the same time. In Fig. a green and red dots indicate positions of corrector and 
BPM to be installed, respectively. At a first (second) BPM, a horizontal and vertical rms 
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FIG. 7: Orbit distortion induced by errors listed in Table |TT] with 100 random seeds. Green and 
red dots mean the position of correctors and BPMs, respectively. 


TABLE II: Errors for orbit correction 


- 

Mialignment [cm] 

Tilt [mrad] 

Field amplitude [%] 

Solenoid 

150 

5.0 

0.05 

Bending 

150 

5.0 

0.05 

Quadrupole 

150 

5.0 

0.05 


orbit sizes are 0.5457 (0.0154) and 0.0113 (0.01) cm, respectively. 

After correcting distorted orbits, transverse orbit trajectories are shown in Fig[^ After 
correction, the rms orbit sizes at first BPM is 0.0190 cm in horizontal plane and 0.0021 cm 
in vertical plane, respectively. At second BPM, those are 0.0027 and 0.0019 cm, respectively. 
The orbit sizes decrease drastically after correction, espically in horizontal plane. During 
correction of distorted orbits for 100 random seeds, average horizontal and vertical kick angle 
of first corrector is 2.51 and 0.25 mrad, respectively. The kick angle of second corrector is 
0.55 and 0.39 mrad which are less than first one. The calculated kick angle of each corrector 
is much less than the mechnical limitation, about 7 mrad. 

The orbit correction is carried out with less number of correctors as shown in Fig. 
Without the second corrector, the rms horizontal and vertical orbit sizes at first (second) 
BPM are 0.0187 (0.0187) and 0.0047 (0.0049), respectively. These values are larger than the 
above results of 2 correctors, however beam sizes at first BPM are much less than the results 
of before correction. The result of orbit correction with different number of correctors is 
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FIG. 8: Horizontal (top) and vertical (bottom) beam centroids after orbit correction with 2 cor¬ 
rectors and 2 BPMs. Green and red dots mean the position of correctors and BPMs, respectively. 
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FIG. 9: Horizontal (top) and vertical (bottom) beam centroids after orbit correction with 1 cor¬ 
rector and 2 BPMs. Green and red dots mean the position of corrector and BPMs, respectively. 


summarized in Table |III[ Finding proper number and position of correcors and BPMs to 
correct the effects of errors will be continued. 


IV. CONCLUSIONS 

We presented the new aspects of significant beam dynamics issues in the main LEBT 
of RAON heavy ion accelerator. The design of MHB was modified newly to improve the 
londitudinal acceptance at the RFQ and SCLl. After changing the design of MHB, the 
rms longitudinal emittance decreased about 26 % at the SCLl entrance. Futhermore, twiss 
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TABLE III: Summary of orbit correction 


# of correctors 

< x{y) >rms at BPMl [cm] 

< x{y) >rms at BPM2 [cm] 

0 

0.5457 (0.0113) 

0.0154 (0.01) 

2 

0.019 (0.0021) 

0.0027 (0.0019) 

1 

0.0187 (0.0047) 

0.0187 (0.0049) 


parameters at the RFQ entrance were matched by using 5 electrostatic quadrupoles, accord¬ 
ingly those were matched less than 0.1 % error with respect to design values. In addition, 
error studies were carried out in the main LEBT. Significant increase of transverse beam 
size appeared between MHB and RFQ with magnet errors, and it could be corrected with 2 
correctors and 2 BPMs. The rms beam size at BPMs decreased from about 0.5 cm to 0.02 
cm after orbit correction within maximum corrector kick angle. The research of error effects 
and finding proper number of correctors will be continued. 
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